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Graves' disease (CD) is an autoimmune thyroid dis- 
order that Is inherited as a complex trait. We have 
aenotyped 77 affected sib-palrs with autoimmune 
thyroid disease for eight polymorphic markers span- 
ning the cytotoxic T lymphocyte antigen-4 (C7L4-4) 
region of chromosome 2q31-q33, and for five mark- 
ers spanning the major histocompatibility complex 
iMHC) region of chromosome 6 P 21. Mon-parametrlc 
analysis showed linkage of GD to the CTLA-A region 
with a peak non-parametric linkage (NPL) score of 
3 43 (P = 0.0004) at the marker D2S117. The propor- 
tion of affected full-sibs sharing zero alleles (zO) 
reached a minimum of 0.113 close to D2S117, giving 
a locus-specific X, for this region of 2.2. Families with 
brother-sister slb-palrs showed a peak NPL of 3.46 
[p- 0 0003, X, > 10) at D2S117, compared with 2.00 
IP- 0 02 X, = 1-8) in the families with only affected 
females, suggesting a stronger influence In families 
with affected males. Association between GD and the 
O allele of the Thr17Ala polymorphism within the 
CTLA-4 gene {CTL.A4A/G) was observed using unaf- 
fected sib controls (P = 0.005). Lesser evidence for 
linkage was found at the MHC locus, with a peak NPL 
score of 1.95 (P = 0.O26), between the markers 
D6S273 and TWFrt. We demonstrate that the CTLA-4 
locus (JL = 2.2) and the MHC locus (X, = 1.6) together 
confer -50% of the inherited susceptibility to GD dis- 
ease in our population. 

INTRODUCTION 

Grabs' disease (GD) is a common organ-specific autoimmune 
disorder, whicb is characterized by thyroid hormone over- 



secretion, diffuse fi oiUe and specific «tad <°^««£ 
(termed thyroid-associated otopathy: TAO). GD affec« 
0.4-0.8% of the female population oyer a life ££). and 
has a concordance rate in monozygonc twins of 20^0% com- 
pared with 5-7% in dizygotic wins or female sibs p-*). 
Familial risk studies combined with local population gvj- 
iencc data have been used to estimate that the excess risk of 
GD to a female sib of u GD proband fc) ^ between 10 and 15 
(2,3). In common with other autoimmune disorders 

s likely to have a complex genetic basis. ^^J*"™ 
gene, each contributing in various degrees to the 
cepiibility. Furthermore, auioimmone hypomyroldism (AW. 
which is the other common maturation of thyroid auio- 
2£* occurs with an increased frequency in GD |k*£. 
.ueeestine that both forms of autoimmune thyroid disease 
SB?) cou d share some susceptibiUty'alleles (10.11). S.mi- 

ariy there is an excess prevalence of GD amongst subject, 
S n'typ' 1 mellilus «DDM) and their relatives sug- 

jSnV *at GD and 1DDM could also share suscepubihty 

" an poTulIuon-Led case^ontrol studies have shown a cons- 
ent Nation of GD with ^^^^^^ 
compatibility complex (MHC) haplotypes ^J^M- 
DOBl*i)20J-DQAl*0SOI) in Caucasian populates 
n 10 13-15V however, these MHC allelic vocations may be 
fSent from mo.« of AH (10). Evidence for linkage o GO to 
MHC has also been found in some populauons (10,15-17) bu. 
SaTbeen dilT.cult .o reproduce (18.19). Other cand.dato 
usecptibility loci, including the thyrotropin receptor immu- 
noglobulin heavy chain (Cm), interleukin I ^ e P lor 
S? cell antigen receptor B and the thyroid hormone receptor 
B Bene have been found to be associated with GD in some pop- 
5£j f studied, but not in others (»«U-25^ntly.^ 
dence for varying degrees of genetic linkage have been 
reported between A1TD and markets on the long arms of 
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Toble 1. Phcnolypes of*fTccte*J iib-pairs wiih autoimmune ihycoid di.wajc 





Graces' disease 
only (GD-GD) 


Mixed Graves' and 
AH (GD-AH^ 


All autoimmune 
thyroid disease 


Full sib-pair.; 


66 


6 


72 


Half sib-pain 


5 


0 


5 


Total 


71 


6 


77 



'Families were .elected on the buaia of two aflected GD sibs. GD-AH sib- 
pairs make up addition*! members of the same familicn. 



chromosomes 14 (14q31) and 20 (20qll) and on the X chro- 
mosome (Xq2J) (19,26,27). These studies await confirmation 
in different populations. 

The cytotoxic T lymphocyte antigen-4 (CTLA-4) is a co- 
stimulatory molecule that is a key negative regulator of T cell 
function (28). The CTLA-4 gene lies on chromosome 2q33, 
and recent studies have demonstrated linkage and association 
of IDDM with markers on 2q3 l-q33 (designated 1DDM7 and 
IDDM12) (29-32). Transmission disequilibrium at alleles pf 
two CTLA-4 polymorphisms defines the IDDM 12 region 
(30,31), and the same alleles haye been associated with GD 
and AH in case-control studies (32-36). However, recent 
genetic linkage studies have failed to confirm CTLA-4 as a GD 
susceptibility locus (19,37). In this study, we have taken an 
affected sib-pair approach to examine the CTLA-4 region of 
chromosome 2q3l-q33 and the MHC region of 6p21 for evi- 
dence of linkage to GD. 

RESULTS 

Chromosome 2q31-q33 linkage analysis. 

We examined the cohort of 77 affected AITD sib-pairs (Table I) 
for linkage to eight polymorphic markers over a 30 cM region of 
chromosome 2q3i-q33, which encompassed the IDDM7 and 
CTLA-4 (IDDM 1 2) regions. Non-parametric analysis with the 
GENEHUNTER package showed a broad region of excess allele 
sharing (up to 65%) amongst affected sibs. When all subjects with 
AITD were designated as affected, the peak multipoint non-para- 
metric linkage (NPL) score was 3.43 (P = 0.0004) at the marker- 
D2S117, which is close to the CTLA-4 region (Fig. 1). Designa- 
tion of only GD cases as affected (7 i sib-pairs) showed a peak 
NPL score of 3.07 (P = 0.001; Fig. 1), suggesting homogeneity 
between GD and AH at this locus. The proportion of the 72 full- 
sibs wiih AITD sharing zero alleles (zO) reached a minimum of 
0.113 close to D2S117, suggesting thai the locus-specific X, for 
this region is 11. The 24 families with affected AITD males 
(brother-sister sib-pairs) showed a peak NPL of 3.46 (JP » 0.0003, 
zO < I x 1(H, \ > 10), compared with 2.00 (P = 0.02. zO » 0.132, 
X, = 1.9} in the 40 families with only affected females, suggesting 
that this locus has a much stronger influence in families with 
affected males. An intrafamilial association analysis, using 
unaffected sibs as controls, showed an excess of the G allele at the 
diallelic CTLA4AIC polymorphism (Thrl7Ala) {P = 0.005), and 
of the 112 mobility unit (mu) allele of CTLA4(AT]n (P = 0.02) in 
GD probands (Table 2). The odds ratio (OR) for the G allele at 
CTLA4A1G in GD probands was 2.01 (95% CI 1.13-4.55; 
P = 0.047) compared with unaffected sibs, and for the GG geno- 
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Figure L Multipoint linkage aiwjyu* of eight markers on chromosome 2q3l- 
q33. (A) Percentage Information con rem over the marker map. (B) The NPL 
score obtained by 'scoring iuT suTected subjects with Graves' disease (solid 
line) and AfTD (dashed line) uaing the GENEHUlsfTER package is shown, 
aeiunst the marker roup on Che *-*xi*. The diallelic polymorphism CTLA4A/G 
\svl the nam* map position as CTLA^fATju. The peak NPL score of 3.43 
occurs ar the marker D2S117, 



type (versus AG and AA) was 4,56 (95% a 1.0(5-22.02: 
P » 0.028). The OR for the G allele in affected GD males com- 
pared with affected female sibs was 1.55 (95% CI 0.6&-3.52; 
P = 0.2). 

Analysis of chromosome 2q3J-^33 linkage conditioned for 
the CTLA4A/G genotype 

Families were weighted according to the CTLA4AIG genotype 
of the proband, and non-pammclric linkage analysis was car- 
ried out with AITD subjects designated as affected. The maxi- 
mum evidence for linkage came from the families whose 
proband carried the AG genotype, with a peak NPL score of 
2.76 occurring at D2SJ 17 (Toble 3). Lesser evidence for link- 
age was found in the families with AA- or GG-carrying 
probands, with peak NPL scores of 2. 14 and 1.36, respectively 
(Table 3). Despite the association between GD and the G allele 
of CTLA4A/G, we find no evidence thai the linkage is confined 
to families whose probands carry this allele. 

MHC linkage analysis 

Examination of five microsatellite markers spanning a 23 cM 
region of chromosome 6p2l encompassing the MHC locus 
showed a modest increase in allele sharing (up to 58%) amongst 
the 71 affected GD sib-pairs, with a peak multipoint NPL score 
of 1 .95 (P * 0.026) occurring at the markers D6S27S and TNFa 
(Fig. 2). Scoring all subjects with AfTD as being affected led to 
a decrease in the peak NPL score ( 1 .40), in keeping with GD and 
AH having different MHC susceptibility alleles. The minimum 
zO value for the 66 full sib-pairs with GD was 0.161 at the 



Tnblc 2. Puir-w.sc wiciadon bciween affected (GD) and un.ffecied sib 
controls at linked chromosome 2 H 1 l-q33 markers 



MurVcr 



Allele 



Graves' 



Unaffected sibs 4 P-«alue* 



U2SSS9 203 mu 


10/12 


5/17 


NS 


P25/77 206 mu 


8/10 


3/15 


NS 


CTIA4AIG C 


20/6 


10/16 


0.005 




14/4 


7/M 


0.020 


l 50mU 


9/7 


3/13 


NS 


■Number of occurence* of candidate aiieles/non-occuirencca (45). 


Wiilwr'a exact teat, corrected tot multiple allelic comparison*. 


Table 3. Analysis of chromosome 2q3 1-^33 Uftkags. condtrioncd for the 


CTLA4AIG genotype of the proband 






Gcnorypw of probands peak NFL score 


/•-value 


No. of families 


CG 




0.087 


13 


AG 2.76 




0.003 


30 


AA 2.14 




0.016 


21 


OG or AG 3.05 




0.001 


43 


AG or AA 3.28 




0.0006 


51 


AH families 3.43 




0.0004 


64 



marker 7WFa, suggesting that the locos-specific X, for the MHC 
region is 1.6. 

Analysis of CTIA-4 linkage conditioned for MHC 
haplotype sharing 

Families were subdivided according to identity by descent 
(JBD) haplotype sharing at the TNFa locus, and the linkage ^to 
the chromosome 2 markers was reexamined. The peak .NFL 
score for sib-pairs sharing t-o MHC haplotypes was 2.38 « 
0.008), compared with 2.19 {P = 0.015) for sibs sharing 1 or 0 
MHC haplotypes. 

DISCUSSION 

We have demonstrated, for the first time, unequivocal cvi- 
dence (P - 0.0004) for linkage of GD to the D2SI 17 region of 
2q33 (Fifi. 1), with a locus -specific \ of 2.2. The type I dia- 
ries \JlDDM7 and the CTLA-4 gene (IDDMJ2) are close to 
the linked region, and it is possible that both may contain sus- 
ceptibility polymorphisms that contribute to the linkage we 
observe Our association analysis, "hich sho^s evidence for 
ititrafamilial allelic association at the 
(CTLA4A/G and CTLA4[AT)«) but not at other loci within this 
region (Table 2). would suggest that a susceptibility polymor- 
phism(s) lies aU or close to. the CTLA-4 locus However, sub- 
group analysis shows that the evidence for linkage to this 
recion i* not confined to families with G allele-carrying 
CTLA4A/C genotypes, such that the effect of this polymor- 
phism alone is not sufficient to explain the observed linkage. 
The strength of the linkage that we observed in the 24 families 
with affected male GD members (P = 0.0003. \ > 10), also 
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l^wn «»in« the mwfcer ord« on the *axU. Tbe peak NPL score of 1.9S 
occurs it the muter* D6S27i and TNFa. 

suggests that this polymorphism(s) may hf e a critical role in 
the susceptibility of males to GD. 

Previous studies of the CTLA4[AT]n marker have not shown 
linkage to Arm; however, both these studies have employed a 
Metric memod of linkage analysis which may be aeo*uveto 
^-specification of background allele frequencies aid the mode 
of inheritance (1 9,37). Funhermore, these studies combined fam- 
ilies from a mixed ethnic background and used a more phenoiypi- 
callv varied population, with many femmes having only un- 
affected subjects (19,37). Thus, our current "^h»™P«J* 
methodological differences from (he CTLA-4 ^ 
reoorted previously. Despite this, genuine differences >n the con- 

between the UK and other populations, as has beer 
demonstrated for the CTLA-4 and other loci, in different IDDM 
populations (6-930.31)- , . _ 

in our population, MHC appears to have ^ >nflu«n« «• 
' GD susceptibility than the CTLA-4 region wuh weaker eW- 
dence of linkage (NPL score 1.95. P « 0.026 «nd • 
4n «:ific X. of 1 6 (Fig. 2). Th s contrasts with the findings m 
S£ MHC has I consistently strong effect in all pop- 
ulations (6-9). However, our finding is not «"«P e ? l ~- " 
>95% of Caucasian IDDM subjects have a 'suscqphble IDDM 
MHC haplotype {DR3IDR4), compared with -50% of conuol 
ISS^JSU only 5<W50% of GD subjectscarry the GD 
susceptible haplotype (DR3). compared with 20-30% of eon 
trol subjects. Thus, MHC has a stronger influence on the devc-l 
opment of IDDM than on GD. Our results also show that 
SSL lmk. 8 c studies of mUed AlTD populatt on, .(mixed 
GD/AH families) (18,19) may have been unable 
age of MHC to ArrD due to allelic heterogeneity between GD 
and AH at this locus (10). . 

Using a multiplicati-c model, we can estimate that the 
CTui Joc^s CK - 2.2) confers 29-34% of the total genetic 
StBCepdbilUy toGD in our population (38). It b thetefore 
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unlUcely ihat any other locus will have an effeci on CD suscep- 
tibility thai is stronger ihan that of CTLA-4. Our data also allow 
us to estimate, for the first time, thai the MHC region (\ = 1 .6) 
confers 17-20% of the genetic susceptibility to GD. Thus, 
taken together, these two loci account for -50% of the inher- 
ited predisposition to GD in our population. 

MATERIALS AND METHODS 

Patients 

Sixty-four families with two or more sibs affected with GD 
(including J 46 with GD, 20 with AH and 12 unaffected subjects) 
were recruited from the north of England and the Lothian region 
of Scotland. GD was confirmed by the finding of biochemical 
hyperthyroidism, with evidence of one of the following: (i) 
significant TAO (American Thyroid Association Class 3 or 
worse) (39); (ii) diffuse increase in thyroid uptake on radio- 
nuclide scan; and (iii) positive serum myrotropin-binding 
inhibitory immunoglobulin antibodies. The cohort of families 
comprised 53 fulUibs both with GD, five GD sib-trios, one GD 
quartet, five half-sibs with GD and six mixed GD/aH fuU-sibs 
(Table 1). Parents (n = 49) and unaffected sibs (n ■ 36) w er e 
studied wherever available. Additional second degree relatives 
had GD (n = 4), AH (n = 3) or were unaffected (n = 5). There 
were 122 female and 24 male GD patients, with a mean age at 
onset of 35.5 years (range 9-67 years). Fifty-four (37%) of the 
GD patients had significant TAO and two had thyroid dcrmc- 
pathy. All members of these sibships were Caucasian, and >95% 
of grandparents were of mainland UK or Irish origin. DNA from 
normal control subjects without evidence or family history of 
autoimmune disease were also obtained from the local popula- 
tion. All studies were carried out with the approval of the 
regional and district ethics committees. 

Genotyping 

The microsatdlite markers were genotyped using fluorcscently 
labelled PCR and resolved on a semi-automated 373 sequencer 
(Applied Btosystems, Foster City ? CA). The primers, except for 
CTLA4(AT}n l were taken from the Genethon genetic linkage map 
(http://ww^.gcnethon.rT/gencthon_.cn.html ). The CTLA4[AT]n 
primer's were: 5-GCC AGT GAT GCT AAA GGT TG-3' and 
5-ACA CAA AAA CAT ACG TGG CTC-3'. Using these prim- 
ers, the 112 mu allele of CTLA4[AT]n is equivalent to the 106 
mu allele described previously (33). The marker map was 
derived from Genethon and the Southampton University data- 
base (htt^://cec^.gcnedcs.sotorLac.uk/public_huTiI/gmap.htmj ). 
The CTLA-4 A/C polymorphism in exon 1 of the CTLA-4 gene 
was amplified using the following primers 5'-CCA CGG CTT 
CCT TTC TCG TA-3' and 5-aGT CTC ACT CAC CTT TGC 
AG-3' followed by digesdon with the restriction enzyme Bstlll 
(Prornega, Southampton, UK) (32). The BsilU digestion assay 
for the CTLA4A/G polymorphism was validated by direct DNA 
sequencing in fWe subjects, as described previously (40). . 

Statistical analysis 

Two-point and multipoint NPL scores and marker information 
content were calculated using the 'score all* function of the 
GENEHUNTER package (41). The minimum proportion of 



full sib-pairs sharing fccro alleles (zO) was calculated for each 
region using MaPMAKER/SIBS (42). Data from the chromo- 
some 2 markers were weighted (0 or 1) for CTLA4AIG allele 
status using the modified GENEHUNTER-plus version 2 soft- 
ware, as described previously (43,44). Tbe population allele 
frequencies for each marker were derived from local Cauca- 
sian controls. Family-based association analysis was per-, 
formed on probands using an unaffected sib as an intrafamiljal 
control (45), and were analysed using Fisher's exact test, re- 
values were Bonferroni corrected for multiple allelic compari- 
sons, except for the markers CTLA4AJG and CTLA4[AT)n, 
where candidate alleles were known. ORs were calculated by 
Woolfs method (46). 

ABBREVIATIONS 

AH, autoimmune hypothyroidism; Aim autoimmune thyroid 
disease; Ct confidence interval; CTLA-4, cytotoxic T lympho- 
cyte antigcn-4; GD, Graves 1 disease; HLA, human leukocyte 
antigen; IDDM, type 1 diabetes raellitus; MHC. major histo- 
compatibility complex; mu, mobility units; NPL, non^parametric 
linkage; OR, odds ratio; TAO, thyroid-associated orbitopathy. 
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